7-2006 



1 



ADMISSIBLE WAVEFRONT SPEEDS FOR A SINGLE SPECIES 
REACTION-DIFFUSION EQUATION WITH DELAY 



Elena Trofimchuk 

Department of Mathematics II, 
National Technical University 'KPI', Kyiv, Ukraine 

Sergei Trofimchuk 

Instituto de Matematica y Ffsica 
Universidad de Talca, Casilla 747, Talca, Chile 

Abstract. We consider equation ut{t, x) = Au{t, x)—u{t, x)+g{u{t—h, x)) (*), when 
g : R+ — > R+ has exactly two fixed points: xi = and X2 = k > 0. Assuming that 
g is unimodal and has negative Schwarzian, we indicate explicitly a closed interval 
C = C{h, g'(0), = [c*, c*] such that {*) has at least one (possibly, nonmonotone) 

travelling front propagating at velocity c for every c £ C. Here c« > is finite and 
c* G R+ U {+oo}. Every time when C is not empty, the minimal bound c« is sharp 
so that there are not wavefronts moving with speed c < c« . In contrast to reported 
results, the interval C can be compact, and we conjecture that some of equations 
(*) can indeed have an upper bound for propagation speeds of travelling fronts. As 
particular cases, Eq. (*) includes the diffusive Nicholson's blowflies equation and the 
Mackey-Glass equation with nonmonotone nonlinearity. 

1. Introduction. In this paper, we study the existence of positive and generally 
nonmonotone travelling waves for a family of delayed reaction-diffusion equations 

Ut{t, x) = Au{t, x) - u{t, x) + g{u{t ~ h, a;)), u{t, x) > 0, x e K™, (1.1) 

which has exactly two non-negative equilibria mi = 0, U2 = k > (so that g{K) = 
K, 5(0) = 0). The nonlinearity g is called the birth function. It is therefore non- 
negative, and generally nonmonotone and bounded. In particular, with g{u) = 
pue'"^ in (|l.l(l . we get the diffusive Nicholson's blowflies equation proposed in |25j . 
Since the biological interpretation of u is the size of an adult population, we will 
consider only nonnegative solutions for l|l.l|l . 

Recently, the existence of travelling fronts connecting the trivial and positive 
steady states in (ll.l|) (and its non-local multi-dimensional generalizations) was 
studied in [3 [5, .20^ 21i I26L |2] by means of two essentially different approaches 
(see also [3] |H1 El for other methods which can be applied to analyze this 
problem) . The first method, which was proposed in (271 131L 133 ' and then further 
substantially developed in I15L I20L 1211 126| . uses positivity and monotonicity 
properties of the integral operator 

{Ax){t) = -^-^ 1 1* e'^'-^^g{x{s - h))ds + e^^'~'^g{x{s - h))ds^ , 
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where A < < /i satisfy ez^ — z — 1 = and e^^/^ = c > is the front velocity. 
As it can be easily observed, the profiles x S C(R, M+) of travelling waves are 
completely determined by the integral equation Ax — x and the first approach 
consists in the use of an appropriate fixed point theorem to A : K ^ K, where 
convex and closed set C {a; : < < x{t) < C C(M,R+) belongs to 

an adequate Banach space (C(IR,K.+ ), | • |). Hence, the existence of travelling fronts 
can be established if we will be able to realize Ma-Wu-Zou reduction which consists 
in choosing profile's subset K of an appropriate Banach space (C(R, IR+),| • |). 
This space should be 'nice' enough to assure the compactness (or monotonicity) 
of A while ip^ have to satisfy A^)^ > <j)^ and A<j)^ < 0+. All these requirements 
are not easy to satisfy, and this why only relatively narrow subclasses of 
(e.g. monotone or quasi-monotone in the sense of I31|) were considered within 
this approach. Here, we would like to mention the recent works I21L 1^ 

which open new interesting possibilities of this method. It should be observed 
that a correct choice of (f>^ is related to the calculation of the minimal speed of 
wavefront's propagation, e.g. see [TTIOTII^I??] 

The second method was proposed in It essentially relies on the fact that, 
in a 'good' Banach space, the Frechet derivative of Mm^^o A along a heteroclinic 
solution -0 of the limit delay differential equation x'{t) = —x{t) + g{x{t — h)) is 
a surjective Fredholm operator. In consequence, the Lyapunov- Schmidt reduction 
can be used to prove the existence of a smooth family of travelling fronts in some 
neighborhood oftp. As it was shown in 5 for the case of this family contains 

positive solutions as well. Furthermore, the method of |31 EI allows to make a 
conclusion about the uniqueness (up to translations) of the positive solution to Eq. 
for every sufficiently small e. Additionally, the following condition established 

in ISJ 

e-''>-rin^^, r'^^g'i^), (1.2) 

and which is sufficient for the existence of positive travelling fronts in for 
g'{0) > 1 and small e > 0, is also 'very close' to the sharp one (see 5, 19 for 
more details and references). However, this analytical approach has one important 
limitation: it can be used only to prove the existence of travelling fronts for suffi- 
ciently small e > 0, or, in other words, for sufficiently large speeds c > 0. On the 
other hand, calculation of the minimal wave speed (and the asymptotic speed of 
propagation) is very important in applications, see discussion about the marginal 
stability conjecture and related questions in [11 1^1^150115^ . 

This work is inspired by the above presented approaches and, especially, by 
[3 |S1 I12L I2UI 1271 I31| . For a broad family of nonlinearities g, we prove that 
Eq. I|l.l|l has positive travelling wave fronts u{t, x) = 4>{ct + v ■ x) provided that 
^'(0) > 1, and c exceeds the minimal velocity: c > c^,{h, g'{0)), and 

^('^'C) = — A)r^^>Y^- (1-3) 

These fronts generally are not monotone: in fact, they can oscillate infinitely about 
the positive steady state. We also show that, for large negative values of s, the wave 
profile (f){s, c) is asymptotically equivalent to an increasing exponential function. 
The minimal speed c^{h,g'{0)) is determined as the unique positive number for 
which the characteristic equation (zjcf' — z — 1 + g'{0) exp{—zh) = has a real 
multiple root. This value of c^,{h, g'{0)) is sharp: there are no travelling waves 
moving at speed c < c*(/i, ^'(0)). 
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Now, notice that < e < £,{h,c) < 1 and £,{h,+oo) = e ^. Thus incquahty 
(|1.3() can be replaced by a stronger one: 

We show that, under the conditions g'{0) > 1 and H1.4|l . the speed interval C is 
infinite, in fact, C — [c*, +00). The latter result is an important supplement to the 
existence conditions proved in [S] since (|1.4|l can viewed as a linearized version of 
()1.2|) where ln(l + w) is replaced with w > ln(l + w), w > 0. 

Next, if is fulfilled and 11.4|l does not hold, we still can prove the existence 
of a compact interval C = [c*,c*] of admissible speeds. Here c* = c*{g'(K),h) is 
determined from the equation 

fi{c*)-x{c*) ^ r^ + r 

In order to state the main result of this work, we need to present our basic hypoth- 
esis: 

(H): Let g e C(M-)_, M+) have only one point of local extremum xm (maximum) 
and assume that g{0) = 0, g{x) > if a; > 0, and 

1 < g'_{0+) = hminf g{x)/x < limsup g{x)/x — .g+(0+) < +00. 

We suppose further that the equation g{x) = x has exactly two fixed points 
and K > 0. Set C2 = g{xM) and let (i < min{g^{xM), xm,k} be such 
that g(Ci) — T^^T^sG[CiX2] di^)- Additionally, g e C^[Ci,C2] and the Schwarz 
derivative Sg is negative for all a; G [Cii C2], x ^ xm- 
{Sg){x) = g"'{x){g'{x))-' - (3/2) (g"(x)(ff'(x))-i)' < 0. 
In the sequel will denote the unique positive number for which the characteristic 
equation {zjcf — z — 1 + g'^{0+) exp{—zh) — has a real multiple root. Also, we 
set C+ = [c;|!",c*], where c* = 00 if H1.4II holds, and c* is determined from (|1.5|) if 
(|1.4I) does not hold and (|1.3|) is true for some positive c. Now we are ready to state 
our main result: 

Theorem 1.1. Assume (H). Then, for every c E C+ = [c+, c*] and v E M™, = 

1, equation has a positive wavefront u(t, x) = (f>(ct + ly ■ x,c). Moreover, 



and cf = {g\.iO+), h) decreases monotonically to as h +00 for every fixed 
(7^(0+) > 1. // additionally g is C^-smooth in some neighborhood of x = 0, then 
the value of — is sharp (minimal speed of travelling waves). Finally, the 
travelling profile (j){t) oscillates about k on every interval [r, +00) provided that the 
equation 

{z/cf ~ z-l + g'{K)exp{-zh)^0 (1.6) 
does not have any root in (—00, 0) U iM. 

We notice again that, in Theorem 11.11 the interval C can be finite exhibiting 
a situation which is not possible for the reaction-diffusion monostable equations 
without delay |7j Theorem 8.3 (ii)] In fact, we conjecture that, for the smooth 
nonlinearity g satisfying (H), it holds C = [c*,c*pt] where c*pt = c*p^{g'[K),h) E 
K U {00} is defined as the largest value for which characteristic equation H1.6|l has 
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exactly one root in the open half plane {5ftz > 0}. In Proposition 15.61 we give 
an argument in favor of this conjecture, establishing that [c,,c*] C [c*,c*p(]. The 
validity of the upper bound c*p( will be discussed in a forthcoming paper. Moreover, 
in view of the Smith global stability conjecture and Lemma l^Tl with Remark 

O below, we believe that C = [c*,+oo) if and only if the scalar delay differential 
equation 

u'{t) = ~u{t) + g{u{t - h)), M > 0, (1.7) 

is globally stable. 

The use of the Schwarz derivative in (H) is motivated by the following proposi- 
tion due to Singer, see I12|ll7j : 

Proposition 1.2. Assume that a : [C*, C*] ~^ [C*i C*]i f G C'^iO'^ b], is either strictly 
decreasing function or it has only one critical point xm ( maximum ) in [C* , C*] • U ihe 
unique fixed point k € [C* , C*] of a is locally asymptotically stable and the Schwarzian 
derivative satisfies {Sa){x) < for all x ^ xm then k is globally asymptotically 
stable. 

It is easy to see that Theorem 11.11 applies to both the Nicholson's blowflies 
equation and the Mackey- Glass equation with nonmonotone nonlinearity |18| . On 
the other hand, the assumption about the negativity of Sg (which requires 
smoothness of g) can be considerably relaxed with the use of a generalized Yorke 
condition introduced in ^ ^1 . 

The structure of this paper is as follows: in the next section, we study the 
dependence of roots to the characteristic equation (|1.6|) in the positive steady state 
of (|l.l|l on the parameter e = c^^. The same for characteristic equation 

(z/c)2-z-l+g'(0)exp(-z/i) =0 (1.8) 

in the trivial equilibrium. R,emark l2.4l of this section contains an estimation of the 
minimal speed c* and shows the monotonicity of c* with respect to h. In Section 3, 
we realize an appropriate Ma-Wu-Zou reduction, indicating a 'good' Banach space 
(C(K,R+), I'D and its subset K. The novelty of approach presented there consists 
in the very simple form of the both bounds (j)^ for K, we believe that our idea can 
simplify selection of upper and lower solutions in other similar situations. Section 
3 contains our first results about the existence of positive bounded solutions to the 
integral equation Ax = x. These solutions are prototypes of wavefronts (in fact, 
following the terminology of |7], they are semi- wavefronts to 0). Before proving 
in Section 5, Theorem 15. II that they are indeed asymptotically constant at +oo, 
we firstly establish their persistence in Section 4. Furthermore, in Remark 15 . 31 and 
Theorem 15 .41 we prove that c* is indeed the minimal travelling speed. Additionally, 
in Section 5 we are concerned with the asymptotic formulae for travelling profile 
(p{t) as t —> —oo (Remark 15.5(1 and discuss the problem of non-monotonicity of ^(t) 
as i — > -foo f Remark 15. 2|l . Finally, in Proposition 15 . 61 we show that our main result 
is consistent with the above stated conjecture about the exact speed interval for 
Eq. jni). 

2. Characteristic equations. In this section, we are concerned with the estima- 
tion of the minimal speed of propagation of wavefronts. For our model, this speed 
is completely determined by the linearization of (|1.1|) about the trivial steady state. 
The other question of interest is the behavior of roots to the characteristic equations 
((1.8(1 and 11.6(1 as functions of parameters e, h. Besides of the use of these results in 
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Theorem ll.il we need them to give a partial justification of our conjectures stated 
in Introduction. Everywhere below, we set 

(f){z) = —z — 1 + a exp(— z/i), "0(2:, e) = ez^ — z — 1 + a exp(— z/i). 

Lemma 2.1. Take a < and h > and assume that JRAj < for every solution Xj 
of the characteristic equation 4>{X) — 0. Then, for every fixed eg > the equation 
"(/"(A, eo) = has only one root Ao in the right semi-plane 5Rz > 0. Moreover, 
Aq e (0, +CX)). 

Proof. The existence of the unique positive root Aq to 7/'(A, eq) — (with some fixed 
eo > 0) is evident. Next, every multiple root zq to this equation has to satisfy the 
system 

eoZg ~ zq ~ 1 + aexp(— zq/i) = 0, 2eoZo — 1 — ahexp{~-zoh) ~ 0, 
which implies 

(coZq - zq - l)h + 2eoZo -1 = 0, a exp{-zoh) = ^ ^ ^° ■ (2.1) 

2 + Aizo 

The first equation of (|2.1() implies that zq is real while the second equation of (|2.1() 
says that zq < 0. Now, it is easy to see that zq is at most of the multiplicity 2, 
and that (zo,eo) is a bifurcation point where two real roots merge and disappear 
as e eo+- 

Hence, each root Xj ^ R of ?A(z, eo) = determines a unique smooth function 
Aj(-) : {aj,(3j) — > C defined on some maximal open interval (aj,f3j) C [0, +oo) 
containing cq and such that Aj(eo) = = 0. Additionally, we claim 

that the arc Aj(-) : {aj,Pj) C can not intersect the imaginary axis from left to 
right. Indeed, we have that 

z^ 

' 2ez-l + /i(ez2-z-l)' 
so that, at the moment of the eventual intersection we have Aj(e*) = iuj and 

nX'jie,) = -u;\l + h + ehu;^)/{{l + h[eu^ + 1))^ + uj^[2e - hf) < 0, (2.2) 
a contradiction. 

Now, if Xj ^ M and Xj{e) were bounded as e — > aj+ > 0, then Xj{aj+) = zq 
would be the negative real root of the multiplicity 2. In consequence, Xj{e) takes 
only real negative values for e > aj , a, contradiction. Therefore, if aj > and 
Xj ^ M, then lime^a^.+ 5R(Aj(e)) = — oo. 

Next, in the event that aj = and {Aj(e„)} is bounded for some e„ 0+, we 
can assume that limAj(e„) exists and satisfies 0(A) — 0. Therefore 5ftAj(e„) < 
implying '^Xj{e) < for all e e {aj,f3j). If aj — 0, and Xj{e) is unbounded as 
e ^ 0+, and 5R(Aj(e)) > in some right neighborhood of 0, then eXj{e) ^ 1 as 
e 0+ so that di.Xj{e) +oo when e ^ 0+. It should be noticed now that, for 
small e > and sufficiently large C > 0, the function iIj{X, eq) has, in virtue of 
the Rouche theorem, only one positive root in the half plane SRz > C. Hence, we 
can conclude that the unique branch Xj{e) with infinite Aj(0+) corresponds to the 
positive real root of 0(A, eg) =0. □ 

Remark 2.2. Analyzing proof of Lemma f2. II we see that if, for fixed a < 0, ft. > 
and for all e > 0, equation ipiX, eo) = has a unique root in the half plane {SRz > 0}, 
then the quasipolinomial 0(A) must be stable. Furthermore, (|2.2|l implies that, if 
for fixed a < 0,/i, eg > 0, the equation 0(A, ep) = has a unique root in the half 
plane {5Rz > 0}, then this property will be maintained for all e > eo. 
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Next proposition, which can be considered as an analog of Lemma l^Tl for positive 
a > 1, is crucial for the calculation of the minimal speed of propagation. 

Lemma 2.3. Assume that a > I. Then there exists a positive sq = eo{a,h) such 
that, for every fixed e e (0,eo), equation ■0(A, e) = has exactly two real roots 
< Ai(e) < A2(e). If e > eq, then ip{X,e) > for all A e M. Moreover, the 
cardinal number N{e) of the set {Xj : ^(A^ , e) = and ^Xj > 0}^^^'' is a decreasing 
function on (0, eo)- Finally, all roots of ip{X,e) = 0, e G (0, Cq) are simple and we 
can enumerate them in such a way that 

• • • < 3fiA4 3fiA3 < Ai < A2. 

Proof. The proof of the existence of such eg = eo(a, ^) is elementary and thus it 
is omitted here. In fact, eo can determined from system (12.111 . where the second 
equation has a unique positive solution zq = zq^o, h) for every a > 1 and ft, > 0. 
Next, for every fixed e > 0, the finiteness of the set {Xj : 5RAj > 0}^^^^ of all roots 
with non-negative real parts to equation il){X, e) = is a well known fact (e.g. see 
P.18]). The monotonicity of N{e) can be deduced from (|2.2|l by repeating the 
same argument as in the proof of Lemma 12.11 The proof of simplicity of roots 
when e G (0,eo) is straightforward (e.g. see Lemma 13 in [S]). Finally, if e > is 
sufficiently small then the last assertion of Lemma [2.31 follows from Lemmas 7,13 
of IH]. Take now the curves Aj(e) and suppose for a moment that ei G (0,eo) is the 
minimal value for which 3fiA4(ei) — 5ftA3(ei) = Ai(ei) < A2(ei). Then the vertical 
line — Ai(ei) contains three simple roots of eiz^ — z — 1 + aexp(— z/i) = 0. On 
the other hand, 3fi(eiz^ — z — 1 + aexp(— z/i)) = with 5Rz = Ai implies that 

ae-^''' cos(/i3z) = e{^z)^ + (Ai + 1 - eXj) = e{^z)^ + ae^^^'*. 

Therefore, we get 

2e(9z)2 + ae-^i^ sm^{0.5h^z) = 0. 
Thus Qz = and, in consequence, A3(ei) = 5RA3(ei) ~ Ai(ei), a contradiction. □ 

Remark 2.4. Set a = g'{0) > 1 and consider zq = zo{h) as function oih. Then, after 
determining Zq(/i) from the second equation of l|2.1|l . we can use the first equation 
of (|2.1|l to get the following differential equation for eo = eo(ft): 

e',{h) = F{h,eo{h)), (2.3) 

where 



, e(2 + zo(/i)) e(2fte - e + 0.5(/i + V4e2 + 4h^e + K^)) ^ 
l + h + hzo{h) h{he + 0.5{h + V4e2 + 4h^e + h^)) 

Since F{h,e) = e/h - (1 - h)/{h{l + zo{h)h + h)), we find that F{h,e) > e/h 
for h > 1 and F{h,e) < e/h for h e (0, 1]. Taking into account the initial value 
eo(l) = (Ina)^^, and using standard comparison theory for differential equations 
and the relation eo(0+) = 0.25/(a — 1), we obtain that, for every a > 1 and h > 0, 



c,{h,a) = l/Veo(/^,a) < min |2Va - 1, y ^| • (2.4) 

Furthermore, F{h,e) > so that c^{h,a) is strictly decreasing in /i, and we see 
from H2.4|l that c*(-|-oo, a) = for every a > 1. Now, solving numerically the initial 
value problem eo(l) = (lna)~^ for scalar differential equation H2.3|l . we can easily 
plot graphs of c* against h, for every fixed g'{0) > 1. 
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3. An application of Ma-Wu-Zou reduction. Throughout this section, xr- (t) 
stands for the indicator (or characteristic function) of M_. Next, foUowing the 
notations of Introduction and Lemma 12.31 for given e e (0, eo) we will denote 
by Ai = Ai(e) and A2 = A2(e) the positive roots of ^{z,e) = 0, while the roots 
ez^ — z — \ — Q will be denoted by A < < /i. Notice that 

\ < Q < Xi < \2 < ^l. 

Also, everywhere in this section, we require that 

(L): g : M+ R+ is bounded and locally linear in some right 5- neighborhood 
of the origin: g{x) — px, x e [0,(5), with p > 1, and that g[x) < px for all 
x>0. 

Assuming all this, we will prove the existence of non-constant positive bounded 
solutions to the equation 

ex"{t) ~ x'{t) x{t) + gixit-- h)) ^0, teR, (3.1) 

satisfying a;(— 00) = 0. Being bounded, each such solution must satisfy the integral 
equation x = Ax: 

x{t) ^ \ I / e^^*-'^g{x{s-h))ds+ r ef"^*-'^ g{x{s - h))ds\ . (3.2) 



e(A*-A) 

As it was shown by Ma, Wu and Zou O EHl EH EH EEl 123 EH ESI , solving 
can be successfully reduced to the determination of fixed points of the integral 
operator A considered in some closed, bounded, convex and A- invariant subset K 
of an appropriate Banach space {X, || • ||). In this paper, the choice of if C A is 
restricted by the following natural conditions: (a) non-zero constant functions can 
not be elements of X; (b) the functions 0+(t) = (5exp(Ait) and 4>^{t) — (5(e^i* — 
e^^*)XR- {t) G K belong to X; (c) the convergence a;„ ^ x on if is equivalent to the 
uniform convergence Xn =^ xq on compact subsets of R. Thus, with this in mind, 
for some p € (Ai,/i), we set 

X = {xeC(R,R):\\x\\^supe-^''^^\x{s)\+snpe~P'\xis)\<oo}; 

s<0 s>0 

K = {x € X : (t)- [t) = (5(e^i* - e^^*)xR_ {t) < x{t) < Se^'' = (/.+ (t), t e M}. 

We will consider also the following operator L, which can be viewed as the 'lin- 
earization' of A along the trivial steady state: 

{Lx){t) = ^ ( / e^(*-"'a;(s - h)ds + e"'^'-''^ x{s - h)ds 

e(^ — A) U^oo Jt 

Lemma 3.1. We have Lcf)^ = 0+. Next, for every il;{t) = {e^'-* — e'^*)xR_(0 G K 
where v G (Ai, A2], we have {Li}j){t) > 'il'{t), t £ R. 

Proof. It suffices to prove that {Ltp){t) > tp{t) for t < 0. Set p = pe~^{ii — A)~^, 
then we have 

{LiP){t) = P e^(*-'V(s - h)ds + e''(*-"'V(s - h)d.^ > 

>p|y" e^^'-'^^j{s-h)ds + J^^ e^(*-")(e^'*'"''^ -e''*'"''^)rfs| > V'(0- 

□ 

Lemma 3.2. Let assumption (L) hold and e e (0, eo)- Then A{K) C K. 
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Proof. We have Ax < Lx < = 0+ for every a: < Now, if for some 

s we have < 4'~{s) < x{s), then s < so that x{s) < 6e^^'^ < 6 implying 
g{x{s)) = px{s) > p(t)~{s). If 0^(si) = then again g{x{si)) > p<p~{si) — 0. 
Therefore Ax > Lcl)^ > (p^ for every x ^ K. □ 

Lemma 3.3. K is a closed, bounded, convex subset of X and A : K —f K is 
completely continuos. 

Proof. Notice that the convergence of a sequence in K amounts to the uniform 
convergence on compact subsets of M. Since g is a bounded function, we have 

1^^^) - e(M - A) (^-^^ 

for every x ^ K . The statement of this lemma follows now from the Ascoli-Arzela 
theorem combined with the Lebesgue's dominated convergence theorem. □ 

Theorem 3.4. Assume (L) and let e e (0, eq)- Then the integral equation ij^.jj)) 
has a positive bounded solution in K . 

Proof. Due to the above lemmas, we can apply the Schauder's fixed point theorem 
to A: K ^ K. □ 

4. Bounded and uniformly persistent solutions of Eq. ()3.1|> . The following 
assumption (B) will be needed throughout this section: 

(B): g : M+ is continuous and such that, for some < Ci < C2, 

1. .9([Ci,C2]) C [Ci,C2] and5([0,Ci]) c [0,C2]; 

2. min^GlCi^Calffls) = 9{Ci); 

3. g{x) > X for x e (0, Ci] and 1 < 5^(0+) < .9^(0+) < 00; 

4. In [0, C2], the equation 17(2;) — x has exactly two fixed points and k. 
In the sequel, we set C = C{[-h, 0],R) and C+ = C{[-h,0],R+). We will write 
x{t,x'Q,(p) for the solution to the initial value problem a;'(0) = x'q, x{s) — (p{s), 
with {xQ,ip) e R X C. We will also use the customary notation xt & C where 
xt(s) = x{t + s), s e [— /i, 0]. In a standard way, one parametric family of maps 
S'* : R X C R X C, t > 0, 

S\x'q,ip) = {x'{t,x'Q,ip),xt{x'Q,ip)) 

associates with autonomous equation (|3.1() the continuous semi-dynamical system 
S'*. If x{t) is a bounded non negative solution of (|3.1() : < x{t) < Mq, i G K, 
then 1)3. 3|l implies that |a;'(t)| < (e(/i — A))~^ max^g[o.7\/„] In consequence, 

{x'{t),xt) : R ^ M X C is a bounded complete trajectory of 5*. Hence, with every 
non- negative bounded solution x : M ^ R+ of Eq. (|3.1|) we can associate a compact 
invariant non-empty w-limit set io{x) C K x C+. Now, if 

< rn = liminf x{t) < limsupx(t) — M < +00, 

t— »+oo t~t + oo 

then, due to the compactness of i-o{x), it holds 

m — min min (/)(s), M — max max (j){s), 

tj}£PrLj(x) se[-/i,0] ipePri^{x) se[-h,0] 

where Pr : R x C+ C+ is the canonical projection on the second factor. This 
construction has the following important consequence: 

Lemma 4.1. [m, M] C g{[m, M]). 



ADMISSIBLE WAVEFRONT SPEEDS 



9 



Proof. Indeed, if there exist M = maxsgRa;(s) and m — inf^gR a::(s), a straightfor- 
ward estimation of the right hand side of (|3.2|l generates AI < ma^rn<x<M gix). 
As long as the maximum M is not reached, due to compactness of Pruj{x) we 
always can find a solution z{t) of H3.1|) such that z(0) = max<,gRz(s) = M and 
infsgRz(s) > m. Therefore, by the above argument, M < mayim<x<M gi^)- The 
inequality m > minrn,<x<M g{x) can be proved in a similar way. Thus we can 
conclude that [m, AI] C g{[m,M]). □ 

From now on, wc will use the following function ^ : R+ — )■ R+, 

/i — A 



fie — Ae 

Lemma 4.2. If x is a solution of the boundary value problem x{a) — xq, x'{b) = 
for if^) . then 

x{b) = ^{b - a) jxo + ^(^7^ ly^^^'^'^^ ~ e''(''-"))g(x(u - 
Proof. It suffices to consider the variation of constants formula for Eq. %\.\\ : 
x{t) = Ae^* + Se^* + -^^j^^ e^(*-^'5(s)ds + j^^ e^(*-^)5(s)ds| , 

where ^(s) := 5(?/(s - ft,)). □ 

Lemma 4.3. Assume (B) and suppose that supg>o.g(s) < ^2- Let x ^ be a non 
negative solution of Eq. yd.l]) . Then we have the following uniform permanence 
property: 

C,i < liminf a;(t) < limsupa;(t) < C2- 

Proof. We follow the main lines of ^| Corollary 12(3)] and ^| Theorem 3.6(a)]. 
The positivity and the upper bound of x can be easily obtained from the integral 
equation (|3.2|1 . so that we need only to prove the first inequality. We proceed by 
contrary, supposing that Ci > p = liminf x(i). Set A — iTifa<x<Ci[g{x) /x] > 1 and 

>+oo ~ ~ 

let Z > be large enough and £ > be small enough to satisfy A(l — £)(!— ^(^)) > 
1 + e. Now, we observe that x can not decrease monotonically on some semi-axis 
[r, +oo). Indeed, in the latter case, x'{t) < and Ci > x{t) > x{+oo) = for 
sufficiently large t and therefore x{t) < x{t — h) < g{x{t — h)) (see (B3)) implying 

ex"{t) = x'{t) + x{t) - g{x{t - h)) < 0. 

This leads us to a contradiction: x{+oo) = —oo. 

In consequence, we can indicate a local minimum point sq such that x{sq) < 
(1 + e)p and x{s) > p{l — e) for all s E [sq — I — h,so]. Then we have that 
g{x{s)) > A(l — e)p for every s G [so ~ ^ ~ h,so]. Indeed, if x{s) < Ci, we obtain 
glxis)) = x{s){g{x{s))/x{s)) > Ax{s) > Ap{l - e). Now, if x(s) € [Ci,C2] then 
g{x{s)) > 5(Ci) > ACi > Ap{l - e), see (B2). 

Therefore, in view of Lemma 14.21 we have 

x{so) = e(0 l^iso - /) + , ^ ^, / " [e^i'^o-i-u) _ g^(so-i-«))^(^(y _ h))du\ > 



^{l^ - A) Jso-l 

I e(M - A) Jso-i J 
= pA{l — ^{l))(l — £)>(! +e)p, a contradiction. 
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□ 

Lemma 4.4. Assume (B) and suppose that sup^>Qg(s) < (2- Let x ^ be 
a non negative solution of Eq. Hc/.1\) such that rh — miseRx{s) < C,i. Then 
limt^_oo x{t) = 0. 

Proof. Set M = sup^gg a:(s), Lemma [4. II guarantees that 

[m,M] C g([TO,M]). (4.1) 

The assumption (B) makes impossible the case to > 0. In consequence, rh — Q and, 
due to Lemma 14.31 either x{—oo) = or 

= liminf a:(i) < hmsupa;(<) = S. 

However, as we wih show it in the continuation, the second case can not occur. 
Indeed, otherwise it is possible to indicate two sequences of real numbers p„ < (?„ 
converging to —00 such that x has local maximum (minimum) at p„ {qn)i x{qn) < 
x{s) < x{pn) for all s € (Pmln) and lima;(p„) = S, limx(g„) — 0. We notice 
that necessarily lim((7„ — p„) = +00, since in the opposite case an application of 
the compactness argument leads to the following contradiction: the sequence of 
solutions x{t — Qn) contains a subsequence converging to a solution ip G C(R, M) 
of Eq. (|3.2|l verifying ■(/'(O) = and V'(io) = S, for some finite to < 0. Hence, 
x{qn — h)> x{qn) and therefore, for sufficiently large n, we get g{x{q„ — h)) > x{qn) 
contradicting to 

tx"{qn) - x'{qn) - x{qn) + g{x{qn - h)) = 0, 
since x" {qn) > 0, x' {qn) — Q- □ 

Theorem 4.5. Assume (B) and suppose that c > c^{g'+{0+)), then Eq. I\3.1]) has 

at least one positive bounded solution x{t) such that 

= lim x{t) < Ci < liminf x(t) < supa;(s) < ^2- 

Proof. First, let us assume that maxa;>o g(a;) = maxj.g[^j^.^] 17(2;). Take positive 
integer k such that kx > g{x) for all a; > and consider the following sequence 

{fcu, for u £ [0, l/(nfc)]; 
1/n, when u € [l/(nfc),inf (^^""^(l/n)] 
g{u), if w > inf g^^(l/n). 

of continuous functions g,, , all of them satisfying hypothesis [L). Obviously, gn 
converges uniformly to g on R+. Now, for all sufficiently large n, Theorem l3.4l and 
Lemma 14.31 guarantee the existence of a positive continuous function Xn (t) such 
that a;„(— 00) = 0, liminft^+00 Xn{t) > Ci, and 

^n(t) = , ^ ( /* e^(*-^)<?„(a;„(s - h))ds + e^^'~''> g^{x^is - h))ds 



eifi - A) 

Since shifted functions a;„(s + a) satisfy the same integral equation, we can assume 
that a;„(0) = 0.5Ci. 

Now, it is clear that the set {a;„} is pre-compact in the compact open topology of 
C(R, M), so that we can indicate a subsequence a;„^. {t) which converges uniformly on 
compacts to some bounded element x e C(R, M). Since linij^+oo g-nj {xnj {s — h)) = 
g{x{s — h)) for every s g R, we can use the Lebesgue's dominated convergence 
theorem to conclude that x satisfies integral equation 13.2|l . Finally, notice that 
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x{Q) = 0.5Ci and thus a;(— oo) = (by Lemma f4. 41) and liminft^+oo x{t) > Ci (by 
Lemma |4.3|) . 

To complete the proof, we stih have to analyze the case when maxx>Q g{x) > 
max^g[(^j 1^2] 5(2^)- However, this cases can be reduced to the previous one if we 
redefine g{x) as §{(^2) for all x > C21 and then observe that sup(gRa;(s) < (2 for 
every solution obtained in the first part of the demonstration. □ 

Remark 4.6. For Eq. Theorem 14.51 considerablv improves the result given in 

Theorem 1.1 from [21]. Indeed, the method employed in [21] needs essentially that 
limsup„^_i_g(5'(0) — g{u)/u)u~'^ is finite for some v E (0, 1] and that g{u) < g'{0)u. 

5. Heteroclinic solutions of Eq. Everywhere in this section, we will 

assume the hypothesis (H) so that all conditions of the weaker (B) will be satisfied. 
Assume that c > cf and let x{t) be a bounded positive non-constant solution of 
Eq. (|3.1I) which existence was established in Theorem 14.51 Set 



Below, we prove that m — M ~ k, once c is taken from the admissible speed interval 
C+. In consequence, the mentioned x{t) is in fact a "travelling profile" for Eq. 1)1. l|l . 

Theorem 5.1. Assume (H). Then, for every c E C+, equation ys.l]) has a positive 
heteroclinic solution x{t). 

Proof. We have proved already that K,m,M £ [Ci,C2] and that [m,M] C g{[m,M]). 
Hence, the sequence of intervals {g''{[m, M])} increases and 



As it can be easily checked, g{[(^X*]) = [C*jC*] G ^ (i^^ f^-ct, [C*7C*] is the global 
attractor of g : [Ci, C2] — > [Cii C2], see 52] for details). In particular, this means that 
g{C*) — C*- Furthermore, there are two possibilities for g : [C*,C*] ~* [C*iC*]: 
(Fl) g is monotone decreasing on [C*, C*]j ^^^d (F2) g is unimodal on [(^*, C^*]. 

Taking into account these characteristics of g and applying Lemma llTTl we obtain 
that each of the following three relations: k < xm, 01 k < m < M, 01 m < M < k 
implies immediately that m — M — k. Therefore, we will concentrate our attention 
only on the case when m < k < M and k > xm (therefore, g'{K) < 0). 

Now, below we will consider only the 'unimodal' case (F2), analogous and simpler 
'monotone' case (Fl) is being left to the reader. 

Let g^^ : [g'^ixAi), (*] ~* [C*,^m] be the inverse of 5 on [(,,a::j\/] and -0 : [C*:C*] ^ 
[xmj C*] 9 k be the inverse of g on [xm, C*]- 

Since Pruj{x) is a compact invariant set, there exist a trajectory y{t) and a 
real h such that y{h) = M. Since k < M < C*, we obtain from H3.1(l that 
?;(0) < ■0(M) e [t, k], where l = min{.gf ^(M), C*} if gi^{M) ^ and l = C* in the 
opposite case. 

Taking into account the boundary condition y'{h) = 0, setting 17(5) := g{y{s — h)) 
and then using Lemma 14.21 we find that 



m = liminf a;(i) < limsupa;(t) = M. 



KM] c nr=o/(KM]) c nr=o/([Ci,C2]) [C*,C1- 



M 




ah)y{0) + {1 ~ ^h)) mux g{x)<^ij{M) + {l-0 max 5(2;). 

xe[m.M] xe[m,AI] 
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Next, if we set ^ = S,{h) and 6{x) = x — we get 

0{M) < (1 -0 max 5(2;). (5.1) 

xe[m,M] 

Analogously, there exists a solution z{t) such that z{h) — m. Then m £ [C*, k] and, 
since z(0) > m > g^^{ni), (when g^^{m) ^ 0) we have from (|3.1|) that z(0) > ?/;( 
In consequence, 

implying 



TO 



m > S,ip{m) + {1 — £,) min g{x), 

xe[m,M] 



e{m)>il-0 min,, 5(2^)- (5.2) 

3:e[m,Af] 

Next, observe that 

9{[c.,c]) = [c* - fc, c* - ^2:m] 3 [(1 - c)c*, (1 - oc] = (1 - e)5([c*, CD 

and 9 is an increasing function, since ip is decreasing. Evidently, the function fog is 
unimodal if the functions g, f arc unimodal and increasing, respectively. Therefore, 
as 

9{k) ev'(K) - (1 - 0^ - (1 - c)ff(«), 

and the function (1 — (^)g is unimodal, we have that 

a = 9-'oi{l-^)g) : [C„C]^ [C*,C1 

is well defined unimodal map with a unique positive fixed point k. Furthermore, 
as it was proved in 12 , the inequality (S*;;) (a;) < implies {Sa){x) < and, for 
(?'(/«) < 0, the inequality |(t'(k)| < 1 amounts to 

{g'{K)f + 1 

From the estimates (|5.1f) . (|5.2|) we obtain that 

[to,M] C a{[m,M]) C cj'^{[m,M]) C ... C (j^{[m,M]) C ... 

Now, for the unimodal function cr, the last chain of inclusions and the two conditions 
|o''('«)| < 1 and {Sa){x) < are sufficient to deduce m = M ~ k, see Proposition 
11.21 Hence, Theorem 15.11 and the existence part of Theorem 11.11 are completely 
proved. □ 

Remark 5.2 (Nonmonotonicity) . The proof of the final assertion of Theorem 11.11 
about nonmonotonicity is in fact contained in the proof of Lemma 16 of |1] and 
therefore is omitted here. Now, the method applied in requires C^-smoothness 
of 5 at K and the hyperbolicity of Eq. H1.6|) . However, as we show in a forthcoming 
paper, the first condition can be weakened and second one can be removed. In fact, 
it suffices to assume that g is a continuous function which is differentiable at n. This 
observation is also valid for the calculation of the minimal speed of propagation, 
see Theorem 15 .41 below. 

Remark 5.3 (The minimal speed of propagation) . The minimal speed c, for reaction- 
diffusion functional differential equations (|1.1() was already calculated in the pio- 
neering work 2;^ of K. Schaaf. See Theorem 2.7 (i) and Lemma 2.5 of the mentioned 
paper (see also [351 Remark 4.1] for time-delayed reaction-diffusion models different 
from Eq. (|1.1|) '). However, the proof of minimality of c, given in I23j seems to be 
incomplete. Indeed, the proof of Proposition 2.6 in 23 uses, in an essential way, the 
hyperbolicity of Eq. (|1.8() . See the last paragraph on p. 591 and the formula (2.25) 
in |23j . Unfortunately, for some c = c' < c*, Eq. I|1.8(l can be non- hyperbolic (in 
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other words, it can have roots on the imaginary axis) and thus still it is necessary 
to disprove the existence of wavefronts propagating at the 'critical' velocity c' > 0. 
We will do this work below, where invoking Proposition 7.2 from we give a 
complete proof for the minimal speed result. 

Theorem 5.4. If g is C^-smooth in some neighborhood of x = 0, then Eq. 
does not have positive travelling wavefronts propagating at the velocity c < c,,. 

Proof. Suppose the theorem were false. Then there exists a travelling profile x : 
M R_|_, a;(— oo) = 0, moving with positive velocity c' < c*. For the convenience 
of the reader, the following part of the proof will be divided in several steps. 
Step (i) First, we suppose that x{t) — 0(exp(Q:i)), t ^ —oo for some a > (this 
always holds when Eq. H1.8|l with c — c' does not have roots \j on the line = 0. 
See Remark l5 .51 below) . Next, Eq. (|3.2I) implies that, for t < h, 

x{i)>-, / e^(*-")<?(2;(s- /i))ds = — — — / e~'"' gixls - h))ds 

and therefore x{t) can not decay superexponentially at — oo. In view of C^- 
smoothness of 5, solution x{t) satisfies the linear homogeneous equation 

ex"{t) - x'{t) - x{t) + 'i{t)x{t -h) = 0, (5.4) 

where ^{t) — g{x{t — h))/x{t — h) — g'{0) + 0(exp(Q:i)). From |22l Proposition 
7.2] (see also Proposition 2.2]), we conclude that there are b> a, S > and a 
nontrivial eigensolution u{t) of the limiting equation 

eu"{t) - v!(t) - u(t) + g\Q)u(t - /i) = (5.5) 

on the generalized eigenspace associated with the (nonempty) set A of eigenvalues 
with 5Re A = 6, such that x(t) = u(t) + 0(exp((6 + (5)f), t —00. On the other 
hand, since c < c*, we know that there are no real negative eigenvalues of H5.5|l : 
hence 3mA 7^ for all A € A. From Lemma 2.3], we conclude that a;(t) is 
oscillatory, a contradiction. 

Step (ii) Next, if a;(— 00) = and x{t) = 0(exp(at)), t —00 does not hold for 
any a > 0, then there exist a sequence i„ ~^ —00 and a real number d > such that 
x'{tn)/x{tn) < d for every n. Indeed, otherwise for every a > we can find t < 
such that x'{t) > ax{t) for every t < t. This implies that {x{t) exp{—at)y > 
for t < T, and therefore a;(t) exp(— at) < a;(r) exp(— ar), t < t. Thus x{t) = 
0(exp(ai)) at t = —00, a contradiction. Since a;(— 00) — 0, without the loss of 
generality we can suppose that x'{t„) > and x{t) < x(tn) for all t < tn- 

Now, since x{t) satisfies (|5.4|) . we conclude that every yn{t) = x{t+tn) / x{tn) > 
is a solution of 

ey"{t) - y'{t) - y{t) + ^{t + t„)y{t - h) = 0, (5.6) 
and y„(0) = 1 = maXf<o yn{t), < 2/'„(0) < d. A partial integration of (|5.t)|) gives 

y'nit) = V'ni^y'' + - f e^'-'^^yAs) - l{s + tMs - h))ds, (5.7) 
e Jo 

from which we deduce the uniform boundedness of the sequence {y'nit)}: 

\y'n{t)\<d + 2g'{Q), t<0,neN. (5.8) 

Together with < y„(t) < 1, t < 0, inequality (|5.8|l implies the compactness of 
sequence yn{t), t < 0, in the compact open topology of C(M_,M+). Therefore, 
by the Arzela-Ascoli theorem, we can indicate a subsequence yn^ (t) converging 
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uniformly to a continuous function y{t) on bounded subsets of M_ and additionally 
we may assume that limj^oo y'n (0) = Vq exists. In consequence, lim„_^oo "/{t+tn) — 
g'{Q) uniformly on R_. Integrating (|5.7() between t and and then taking the limit 
as j ^ cx) in the obtained expression, we establish that y{t) satisfies Eq. H5.5|l . 
Additionally, y{0) = l,y'(0) ^ y'o € [0, d] and < y{t) < 1, t < 0. Moreover, 
we can prove that y{t) > for all t < 0. Indeed, if y{s) = with s < then 
y'{s) = 0,y"{s) > so that if^ implies y"{s) = and y{s - h) = 0. Setting 
in Lemma 14.21 aix) — g'{Q)x, a — s — h and b = s, we find that y{t) = for all 
t G [s — 2h, s — h]. Thus y{t) = for all t > s — 2h which is not possible because of 

y(o)-i.^ 

Summing up, we have proved that Eq. 1)5.5(1 has a bounded positive solution 
on M_. As it was established in [161 Lemma A.l], this solution does not decay 
superexponentially. In consequence, arguing as in the last paragraphs of Step (i), 
we finalize the proof of Step (ii) and Theorem l5.4l □ 

Remark 5.5 (Asymptotic formulae). Under conditions c > c* as well as of C^- 
smoothness of g and hyperbolicity of H3.1|l at zero, we can repeat Step (i) above to 
derive asymptotic formulae for heteroclinic solutions. Indeed, a;(— oo) = together 
with (|3.2() implies a;'(— oo) = 0. Hence, the solution {x{t),x'{t)) belongs to the 
unstable manifold of the trivial equilibrium of Eq. ((3.1|) . Thus \x{t) \ + — 
0(e"*), t — > — CX3, for some a > (see |21 Section VIII. 4]). On the other hand, 
we know that x{t) is not decaying superexponentially at —oo. Since x{t) satisfies 
the homogeneous equation H5.4|) . due to j23 Proposition 7.2], we get, for some ki 
and every smaU 6 > Q, that x{t) = fcie^^* + 0(e(^i+*)*) = 0(e^i*), t -> -oo. In 
consequence, x{t) satisfies the inhomogeneous equation 

ex"{t) - x'{t) - x{t) + g'{Q)x{t - h) ^ H{t), t e K, 

where H{t) = g'iO)x{t - h) - g{x(t - h)) = 0(exp(2Ait)) as t -oo. Now, from 
[221 Proposition 7.1] and Lemma we deduce that there exist ki,k2, kl + k^^ 0, 
such that, for every fixed (5 > 0, 

r /cie^i* + fee^'* + 0(6(2^1-5)*), t -> -oo, if 2Ai > A2, ki ^ 0; 
x{t)^\ /fcie^i* + 0(6(2^1-5)*), t-> -00, if 2Ai < A2, /fci 7^0; 

\ fee^"* + 0(e(2^2-5)t)^ I _^ _oo, if ki ^0, 0. 

Hence, taking an appropriate and i G {1,2}, we obtain the following asymptotic 
formula for travelling wavefront profiles: 

a;(i - io) = exp(A,t) + 0(exp((2A, - (5)t)), i ^ -00. (5.9) 

To get a similar formula for a;'(i — io)i it suffices to differentiate (|3.2(l and then use 
((5.9() in the right hand side of the obtained integral representation for x'{t — to), 
see ^ Theorem 14] for details. 

Finally, in Proposition 15.61 below we prove that c* < c*p^. In this way we show 
that the existence result of Theorem II. H is consistent with the conjecture about the 
exact speed interval for Eq. 1)1.1(1 that was proposed in Introduction. 

Proposition 5.6. Let us assume i)5.5() and set a — g'{K). Then the characteristic 
equation 1^1. t^] has only one (real) root Aq in the right half plane SRz > 0. 

Proof. On the contrary, let us suppose that Eq. ((1.8(1 has a pair Ai^2 of complex 
conjugate roots with 5RAi,2 > for some fixed values /lo, 5'(k) and ep. Notice that, 
by virtue of Lemma l2. II in this case 

C - ^{ho) > {g'{^)r + 9'i^))/{{9'i^) f + 1) > 



ADMISSIBLE WAVEFRONT SPEEDS 



15 



Here, we use the fact that (|1.4f) imphes ^Xj < for every solution Xj of the 
characteristic equation (p{X) = 0. See |12[I19| . 

Now, keeping eq and g'{K) fixed and considering h G (0, /iq], we define hi < h2 
as the minimal value of h for which equation H1.8|l has roots fi — /i(/ii) — ±iw on 
the imaginary axis. Since 

^dt4^ = u;^{l + 2e + 2eW)/{{l + h{eu^ + 1))^ + u\2e - hf ) > 0, 
ah 

by the Hopf bifurcation theorem 13 , Eq. has a supercritical Hopf bifurcation 
dX h = hi. Therefore, for all h G {hi, ho) close to hi, equation (|3.1I) has a stable 
periodic solution bifurcating from the equilibrium k. By Theoreni ll.il this implies 
immediately that 

a contradiction since ^ — ^(h) is a strictly decreasing function. □ 
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